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Physics beyond the standard model might imply the cosmological production of particles with grand 
unification scale energies. Nucleons and 7-rays from such processes are candidates for the cosmic 
rays observed beyond 100 EeV (10 20 eV). Using a new particle propagation code, we calculate the 
neutrino fluxes predicted by such scenarios if consistency with the observed cosmic ray flux and the 
universal 7-ray background at 1 — 10 GeV is required. Flux levels detectable by proposed km 3 scale 
neutrino observatories are allowed by these constraints. Bounds on or detection of a neutrino flux 
above ~ 1 EeV would allow neutrino astronomy to probe grand unification scale physics. 
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The recent detection of ultra-high energy cosmic 
rays (UHE CRs) with energies above 100 EeV |j|§ has 
prompted an intensive discussion in the literature on the 
nature and origin of these particles ||-|| . At these ener- 
gies, nucleons and heavy nuclei are subject to photopion 
production on the cosmic microwave background (CMB) 
(the Greisen-Zatsepin-Kuzmin effect Q), and to photo- 
disintegration via the giant dipole resonance j7j, respec- 
tively, which limits their range to less than about 30 Mpc. 
Together with the fact that protons above 100 EeV are 
expected to be deflected by only a few degrees in the 
galactic and extragalactic magnetic fields j|, this puts 
severe constraints on models explaining these events by 
acceleration of charged particles in active galactic nuclei 
(AGN) or other known astronomical objects. 

In a more speculative class of so-called "top- 
down" (TD) models, predominantly 7-rays, electrons 
(positrons), and neutrinos are produced directly at UHEs 
by the cascades initiated by the decay of a supermassive 
elementary "X" particle associated with some grand uni- 
fied theory (GUT), rather than being accelerated. The X 
particles are usually thought to be released from topolog- 
ical defects such as cosmic strings, domain walls and mag- 
netic monopoles left over from GUT phase transitions. 
Models of this type predict spectra which are consider- 
ably harder and extend much further beyond 100 EeV 
than shock acceleration spectra and therefore appear as 
a natural alternative to explain the highest energy events 
observed Jq] . The predicted flux level is very model de- 
pendent but, in at least some of the scenarios consid- 
ered [ fTpf , is consistent with the observed UHE CR fluxes. 

Several signatures which might distinguish top-down 



scenarios from more conventional "bottom-up" accel- 
eration models have been suggested in the literature: 
above ~ 100 EeV, domination of the UHE CR flux by 
a hard 7-ray component pd[ | , which potentially produces 
a "gap" |Q in the UHE CR spectrum, and at lower en- 
ergies ~ 10 EeV, a 7-ray to total cosmic ray flux ratio on 
the order of 10% pl| . Another potentially important sig- 
nature of TD models is their prediction of a large cosmo- 
logical neutrino flux || which above 100 EeV can dom- 
inate by far contributions expected from other sources 
such as AGN JlJ] . TD scenarios also predict a substantial 
7-ray flux at lower (GeV) energies that is constrained by 
the observed background levels In this Let- 

ter, we use a new propagation code 18 to establish a re- 
lation between the UHE neutrino flux predictions within 
TD scenarios and the constraints arising from the corre- 
sponding 7-ray flux predictions around 1 GeV and above 
~ 10. EeV. For km 3 scale neutrino observatories, proto- 
types of which are AMANDA, DUMAND, Baikal, and 
NESTOR (l9) , we estimate the UHE neutrino event rates 
predicted by a TD scenario which explains UHE CRs and 
give the maximum event rates consistent with the 7-ray 
flux constraints. We also estimate rates for deeply pen- 
etrating showers for fluorescence detectors such as the 
High Resolution Fly's Eye (HiRes) (2^] and the Telescope 
Array pl[ , and for horizontal neutrino events that might 
be observed in the proposed Pierre Auger project p2| . 

The X particles released from topological defects via 
physical processes such as collapse or annihilation could 
have GUT scale masses mx up to ~ 10 16 GeV. We as- 
sume that the X particles quickly decay into leptons 
and quarks of comparable energy. We take the primary 



1 



lcpton produced in a decay to be an electron with en- 
ergy mj/2. This lepton was not included in prior cal- 
culations but can significantly contribute to the 7-ray 
flux predicted at low energies. The quarks produced 
in the decays interact strongly and fragment into jets 
of hadrons typically containing 10 4 — 10 5 mesons and 
baryons. Given the X particle production rate, dnx/dt, 
the effective production spectrum for particle species 
a (a — 7, N, e , v) via this hadronic channel can be 
written as Q a (E,t) = (dnx /dt)(2/mx)(dN a /dx), where 
x = 2E/mx, and dN a /dx is the relevant effective frag- 
mentation function. For the total hadronic fragmentation 
function dNh/dx we use solutions of the QCD evolution 
equations in modified leading logarithmic approximation 
which provide good fits to accelerator data at LEP ener- 
gies [£3| . We assume that about 3% of the total hadronic 
content consists of nucleons and the rest is produced as 
pions and distributed equally among the three charge 
states. The effective production spectra of 7-rays, elec- 
trons, and neutrinos are then determined from the pion 
decay spectra. The X particle production rate is assumed 
to be spatially uniform and in the matter-dominated era 
can be parametrized as dnx/dt cx t~ i+p , where p 
depends on the specific defect scenario. A network of 
ordinary cosmic strings 0] and annihilation of mag- 
netic monopole-antimonopole pairs fl(]| are represented 
by p — 1, whereas p = 2 corresponds to a constant co- 
moving production rate. As there is considerable uncer- 
tainty in the TD model physics, we will not explicitly con- 
sider any specific TD scenario. Rather, we parametrize a 
TD model by the typical X particle mass, mx, the X par- 
ticle production rate at zero redshift, (dnx / dt)(z = 0), 
and the cosmological evolution of this rate, determined 
by p. 

The shapes of the UHE nucleon and 7-ray spectra pre- 
dicted within TD models are "universal" in the sense 
that they depend only on the physics of X particle de- 
cay. This is because at UHEs nucleons and 7-rays have 
attenuation lengths in the cosmic microwave background 
(CMB) which are small compared to the Hubble scale. 
Cosmological evolutionary effects which depend on the 
specific TD model are therefore negligible. Since the re- 
sulting spectra tend to be harder than any other compo- 
nents from acceleration sources, TD mechanisms could 
contribute to the flux dominantly above ~ 100 EeV, but 
negligibly in the range 10 14 eV — 10 EeV. 

In contrast to the universality of UHE spectral shapes, 
the predicted 7-ray flux below ~ 10 14 eV (the thresh- 
old for pair production of photons on the CMB) and the 
predicted neutrino flux depend on the total energy re- 
lease integrated over redshift and thus on the specific 
TD model. Compared to acceleration scenarios, this en- 
ergy release can be substantial, especially at high red- 
shifts where conventional sources such as galaxies are not 
expected to contribute. The production of UHE parti- 
cles in TD scenarios is therefore subject to a variety of 



constraints of mostly cosmological nature. Electromag- 
netic (EM) energy injected into the universe above the 
pair production threshold on the CMB is recycled into 
a generic cascade spectrum below this threshold on a 
time scale short compared to the Hubble time. This can 
have several potentially observable effects [|l6| such as 
modified light element abundances due to 4 He photodis- 
integration [^|| and distortions of the CMB. Comparison 
with observational data already rules out the class of TD 
scenarios corresponding to p — |l6| ] to which certain 
superconducting cosmic string models belong pG| . 

In addition, for a TD scenario to be viable, its pre- 
dicted spectrum must be consistent with all flux mea- 
surements and limits available for energies between ~ 
100 MeV and a few 100 EeV. Observational data on the 
universal 7-ray background in the 1 — 10 GeV region [^7| , 
to which the generic cascade spectrum would contribute 
directly, turn out to provide an important constraint. 
Since the UHE 7-ray flux is especially sensitive to cer- 
tain astrophysical parameters such as the extragalac- 
tic magnetic field (EGMF), a reliable calculation of the 
predicted spectral shapes requires numerical methods. 
To this end we recently performed extensive numerical 
simulations for the propagation of cxtragalactic nucle- 
ons, 7-rays and electrons (positrons) in the energy range 
10 8 eV < E < 10 25 eV, details of which are given in 
Refs. The calculations take into account all the 

relevant interactions with the (redshift dependent) uni- 
versal low energy photon background in the radio, mi- 
crowave and optical/infrared regime. For nucleons this 
includes production of e" 1 pairs (for protons) and mul- 
tiple pions. The 7-rays interact via production of single 
and double pairs, and electrons via inverse Comp- 
ton scattering and triplet pair production. In addi- 
tion, electrons can suffer significant synchrotron losses 
in EGMFs. We assume a flat universe, a Hubble con- 
stant of Hq = 75 km sec _1 Mpc~ 1 , and zero cosmological 
constant. 

Neutrino fluxes within TD scenarios have been calcu- 
lated before in the literature: Ref. [^9| contains a dis- 
cussion of the (unnormalized) predicted spectral shape 
and Ref. |}(| computes the absolute flux predicted by 
specific processes such as cusp evaporation on ordinary 
cosmic strings. Ref. I| and Ref. |u] also calculate ab- 
solute fluxes within scenarios such as the one discussed 
above, with Ref. |3l| discussing in detail neutrino propa- 
gation and experimental issues related to the HiRes and 
the Telescope Array experiments. However, none of these 
discussions take into account the cosmological constraints 
on TD models, in particular from the low energy 7-ray 
spectrum. 

Fig. [l] shows flux predictions from our code for a typical 
TD scenario for an EGMF of 10 -12 G. In this scenario, 
the UHE CR events are assumed to be cascade 7-rays 
produced by the decay of X particles whose decay rate 
has been normalized accordingly. Although 7-ray pri- 
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maries might be somewhat disfavored [Q, this is consis- 
tent with the currently unknown UHE CR composition. 
The 7-ray and nucleon fluxes shown in Fig. ^ are con- 
sistent with observational estimates of the integral flux 
above 300 EeV [|]J|] , with a likelihood significance above 
50% for E > 100 EeV, and with all data at lower energies. 
Interestingly, the constraints arising from the diffuse 7- 
ray background observed at 1 — 10 GeV are somewhat 
less stringent than earlier analytical estimates |ll],[l(| . As 
a consequence, for EGMFs < 10 -9 G TD scenarios with 
X particle masses as high as I0 16 GeV can still be viable 
models of UHE CRs @. As in Ref. g, the accom- 
panying neutrino fluxes which are also shown in Fig. [l] 
were calculated using the absorption cross section of UHE 
neutrinos in the thermal neutrino background. The re- 
sulting fluxes satisfy present bounds from the Frejus de- 
tector |32| within about five orders of magnitude. With 
fluxes of this level, neutrinos are unlikely UHE CR can- 
didates because of their low interaction probability in the 
atmosphere J33|. 

UHE neutrinos can produce muons in ordinary matter 
via charged current reactions with nucleons (Ns) |54],[35| . 
The most recent calculation of the corresponding cross 
section can roughly be parametrized by <j v n{E) ~ 2.82 x 
10~ 32 cm 2 (£/I0EeV) a402 for E > I PeV where E is the 
neutrino energy |55| . We note that uncertainties of the 
vN cross sections above 10 EeV from extrapolation of 
QCD evolution of up to a factor of 10 translate into cor- 
responding uncertainties in the predicted rates. For an 
(energy dependent) ice or water equivalent acceptance 
A(E) (in units of volume times solid angle), one can ob- 
tain an approximate expected rate of UHE muons pro- 
duced by neutrinos with energy > E, R(E), by multiply- 
ing A(E)a U N(E)ri}i 2 o with the integral muon neutrino 
flux ~ Ej Ufi . Here, iih 2 o is the nucleon density in wa- 
ter. The neutrino energy and arrival direction can be re- 
constructed from the observed muon bremsstrahlung and 
the track geometry. Alternatively, one could use acoustic 
detection methods. The backgrounds are in general ex- 
pected to be small [119 36 1. The rate prediction from the 



model shown in Fig. |1 can be written as 



R{E) ~ 6 x 10 
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around E ~ 10 EeV. Above ~ 100 EeV the corresponding 
fluxes would dominate all present model predictions for 
AGN neutrino fluxes [H as well as the flux of "cosmo- 
genic" neutrinos produced by interactions of UHE CRs 
with the universal photon background j37j]3£||3l|] . Note 
that the p — model in Ref. S which implies neutrino 
event rates about three orders of magnitude higher than 
Eq. (|l|) is ruled out because of overproduction of low en- 
ergy 7-rays [jl3|. We also note that UHE neutrinos can 
produce lepton pairs on the thermal neutrino background 



via the Glashow resonance |39| whose decay products 
in turn contain secondary neutrinos. As was shown in 
Ref. [^9), for mx ^ 10 24 eV this effect, which was not 
taken into account in our simulation, can increase neu- 
trino fluxes around 100 EeV by factors of a few which 
makes our estimate Eq. (^) conservatively low. Our flux 
estimates are further reduced compared to Ref. |u| for 
similar scenarios by our normalization procedure which 
assures consistency of predicted 7-ray and nucleon fluxes 
with observational data at all energies. 
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FIG. 1. Predictions for the differential fluxes of 7-rays 
(solid line), nucleons (dashed line), and (or £> M ), v e (or i> e ) 
(thin solid lines in decreasing order) by the TD model char- 
acterized by p = 1, mx = 2 x 10 25 eV, assuming an EGMF 
of 10 -12 G. Also shown are the combined data from the Fly's 
Eye jj| and the AGASA § experiments above 10 EeV (dots 
with error bars), piecewise power law fits to the observed 
charged CR flux (thick solid line) and experimental upper 
limits on the 7-ray flux at 1 - 10 GeV from EGRET data Q 
(dotted line on left mar gin) . The arrows indicate limits on 
the 7-ray flux from Ref. [Eq] . 

For detectors based on the fluorescence technique such 
as the HiRes p(l and the Telescope Array [^lj, the sen- 
sitivity to UHE neutrinos is often expressed in terms of 
an effective aperture a(E) which is related to A(E) by 
a(E) = A(E)a V N (E)nn 2 o . For the cross sections from 
Ref. |3j|, the apertures given in Ref. |^(| for the HiRes 
correspond to A(E) ~ 3 km 3 x 27rsr for E > 10 EeV 
for muon neutrinos. The expected acceptance of the 
proposed Pierre Auger project for horizontal UHE neu- 
trino induced events is ~ 20 km 3 sr at 10 EeV and ~ 
200km 3 sr at 10 23 cV Q. We conclude that detection 
of neutrino fluxes predicted by scenarios such as the one 
shown in Fig. [I] requires running a detector of acceptance 
> 100 km 3 x 2tt sr over a period of a few years. Again, 
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the backgrounds seem to be negligible |u| . 

A more model independent estimate for the average 
event rate R{E) can be made if the underlying scenario 
is consistent with observational nucleon and 7-ray fluxes 
and the bulk of the energy is released above the pair pro- 
duction threshold on the CMB. Let us assume that the 
ratio of energy injected into the neutrino versus EM chan- 
nel is a constant r. As in Fig. [I], cascading effectively re- 
processes most of the injected EM energy into low energy 
photons whose spectrum peaks at ~ lOGeV pi] ]. Since 
the ratio r remains roughly unchanged during propaga- 
tion, the height of the corresponding peak in the neu- 
trino spectrum should roughly be r times the height 
of the low-energy 7-ray peak, i.e., we have the condi- 
tion max£ [S 2 ji/ M (E)] ~ rmaxg \E 2 j 1 {E)\ . Imposing 
the observational upper limit on the diffuse 7-ray flux 
around lOGeV shown in Fig. ||, max^ \E 2 j Vfi (E)] < 
2 x 10 3 r eVcm _2 sec _1 sr _1 , then bounds the average neu- 
trino rate at all energies E > 1 PeV to be 



R{E) < 0.34 ) 
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For r < 10 5 (10EeV/i?) this bound is consistent with the 
constraint from the Frejus experiment ]32] |. In typical TD 
models such as the one discussed above, r ~ 0.3. How- 
ever, mechanisms with r » 1 could induce appreciable 
event rates above ~ 1 EeV in a km 3 scale detector. A 
detection would thus open the exciting possibility to es- 
tablish an experimental lower limit on r. We stress that 
Eq. ^ holds regardless of whether or not the underlying 
TD mechanism explains the observed UHE CR events. 

The transient event rate could be much higher than 
Eq. (|2|) in the direction to discrete sources which emit 
particles in bursts. Corresponding pulses in the UHE 
nucleon and 7-ray fluxes would only occur for sources 
nearer than ~ 100 Mpc and, in case of protons, would 
be delayed and dispersed by deflection in galactic and 
extragalactic magnetic fields 42 . The recent observa- 



tion of a possible correlation of CRs above ~ 40 EeV 
by the AGASA experiment 0] might suggest sources 
which burst on a time scale tb <C 1 yr. A burst flu- 
ence of ~ r [A(E) /lkm 3 x 27rsr] (E/10 EeV)" - 6 neutri- 
nos within a time tb could then be expected. Associated 
pulses could also be observable in the GeV — TeV 7-ray 
flux if the EGMF is smaller than ~ 10" 15 G in a signifi- 
cant fraction of extragalactic space [44] ]. 

In conclusion, by using a new particle propagation code 
we have given conservative estimates of the UHE neutrino 
flux above 1 EeV which is predicted by a typical TD type 
scenario of UHE CR origin. We demonstrated that the 
constraint imposed by requiring that TD scenarios do not 
overproduce the measured universal 7-ray background at 
1 — 10 GeV implies an upper limit on these neutrino fluxes 
which only depends on the ratio r of energy injected into 



the neutrino versus EM channel, and not on any specific 
TD scenario or even a possible connection to UHE CRs. 
For r ^ 1, neutrino fluxes near this upper limit are poten- 
tially detectable by a km 3 scale neutrino observatory. A 
detection of the UHE neutrino flux might establish an ex- 
perimental lower limit on r and thus allow important in- 
sights into new fundamental physics near the GUT scale. 
Due to the increase of the vN cross section with energy, 
neutrino event rates tend to decrease less strongly with 
energy than UHE CR event rates and the spectral shape 
and cutoff of the neutrino flux (and thus mx) might be 
more easily accessible. A non-detection with more strin- 
gent upper limits would also be useful since it could elim- 
inate large classes of TD models of UHE CR origin. For 
example, failing to detect neutrinos above ~ 10 EeV with 
an exposure A ■ t would rule out scenarios of the type 
shown in Fig. [l] for r > (100 km 3 2n sr yr / A ■ t) . Neutrino 
astronomy might thus be connected to new fundamental 
physics. 
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